Determination of horizon size in state-based peridynamics

Peridynamics is based on integro-differential equations and has a length scale parameter called horizon which gives peridynamics a non-local character. Currently, there are three main peridynamic formulations available in the literature including bond-based peridynamics, ordinary state-based peridynamics and non-ordinary state-based peridynamics. In this study, the optimum horizon size is determined for ordinary state-based peridynamics and non-ordinary state-based peridynamics formulations by using uniform and non-uniform discretisation under dynamic and static conditions. It is shown that the horizon sizes selected as optimum sizes for uniform discretisation can also be used for non-uniform discretisation without introducing significant error to the system. Moreover, a smaller horizon size can be selected for non-ordinary state-based formulation which can yield significant computational advantage. It is also shown that same horizon size can be used for both static and dynamic problems

Positive approximations of the inverse of fractional powers of SPD M-matrices

This study is motivated by the recent development in the fractional calculus and its applications. During last few years, several different techniques are proposed to localize the nonlocal fractional diffusion operator. They are based on transformation of the original problem to a local elliptic or pseudoparabolic problem, or to an integral representation of the solution, thus increasing the dimension of the computational domain. More recently, an alternative approach aimed at reducing the computational complexity was developed. The linear algebraic system $\cal A^\alpha \bf u=\bf f$, $0< \alpha <1$ is considered, where $\cal A$ is a properly normalized (scalded) symmetric and positive definite matrix obtained from finite element or finite difference approximation of second order elliptic problems in $\Omega\subset\mathbb{R}^d$, $d=1,2,3$. The method is based on best uniform rational approximations (BURA) of the function $t^{\beta-\alpha}$ for $0 < t \le 1$ and natural $\beta$. The maximum principles are among the major qualitative properties of linear elliptic operators/PDEs. In many studies and applications, it is important that such properties are preserved by the selected numerical solution method. In this paper we present and analyze the properties of positive approximations of $\cal A^{-\alpha}$ obtained by the BURA technique. Sufficient conditions for positiveness are proven, complemented by sharp error estimates. The theoretical results are supported by representative numerical tests

Integrating Peridynamics with Material Point Method for Elastoplastic Material Modeling

© Springer Nature Switzerland AG 2019. We present a novel integral-based Material Point Method (MPM) using state based peridynamics structure for modeling elastoplastic material and fracture animation. Previous partial derivative based MPM studies face challenges of underlying instability issues of particle distribution and the complexity of modeling discontinuities. To alleviate these problems, we integrate the strain metric in the basic elastic constitutive model by using material point truss structure, which outweighs differential-based methods in both accuracy and stability. To model plasticity, we incorporate our constitutive model with deviatoric flow theory and a simple yield function. It is straightforward to handle the problem of cracking in our hybrid framework. Our method adopts two time integration ways to update crack interface and fracture inner parts, which overcome the unnecessary grid duplication. Our work can create a wide range of material phenomenon including elasticity, plasticity, and fracture. Our framework provides an attractive method for producing elastoplastic materials and fracture with visual realism and high stability

Towards an Efficient Finite Element Method for the Integral Fractional Laplacian on Polygonal Domains

We explore the connection between fractional order partial differential equations in two or more spatial dimensions with boundary integral operators to develop techniques that enable one to efficiently tackle the integral fractional Laplacian. In particular, we develop techniques for the treatment of the dense stiffness matrix including the computation of the entries, the efficient assembly and storage of a sparse approximation and the efficient solution of the resulting equations. The main idea consists of generalising proven techniques for the treatment of boundary integral equations to general fractional orders. Importantly, the approximation does not make any strong assumptions on the shape of the underlying domain and does not rely on any special structure of the matrix that could be exploited by fast transforms. We demonstrate the flexibility and performance of this approach in a couple of two-dimensional numerical examples

A peridynamic based machine learning model for one-dimensional and two-dimensional structures

With the rapid growth of available data and computing resources, using data-driven models is a potential approach in many scientific disciplines and engineering. However, for complex physical phenomena that have limited data, the data-driven models are lacking robustness and fail to provide good predictions. Theory-guided data science is the recent technology that can take advantage of both physics-driven and data-driven models. This study presents a novel peridynamics based machine learning model for one and two-dimensional structures. The linear relationships between the displacement of a material point and displacements of its family members and applied forces are obtained for the machine learning model by using linear regression. The numerical procedure for coupling the peridynamic model and the machine learning model is also provided. The numerical procedure for coupling the peridynamic model and the machine learning model is also provided. The accuracy of the coupled model is verified by considering various examples of a one-dimensional bar and two-dimensional plate. To further demonstrate the capabilities of the coupled model, damage prediction for a plate with a pre-existing crack, a two-dimensional representation of a three-point bending test, and a plate subjected to dynamic load are simulated

Derivation of dual horizon state-based peridynamics formulation based on euler-lagrange equation

The numerical solution of peridynamics equations is usually done by using uniform spatial discretisation. Although implementation of uniform discretisation is straightforward, it can increase computational time significantly for certain problems. Instead, non-uniform discretisation can be utilised and different discretisation sizes can be used at different parts of the solution domain. Moreover, the peridynamic length scale parameter, horizon, can also vary throughout the solution domain. Such a scenario requires extra attention since conservation laws must be satisfied. To deal with these issues, dual-horizon peridynamics was introduced so that both non-uniform discretisation and variable horizon sizes can be utilised. In this study, dual-horizon peridynamics formulation is derived by using Euler–Lagrange equation for state-based peridynamics. Moreover, application of boundary conditions and determination of surface correction factors are also explained. Finally, the current formulation is verified by considering two benchmark problems including plate under tension and vibration of a plate